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Resume 

In this paper we perform an exact study of “Quantum Fidelity” (also 
called Loschmidt Echo) for the time-periodic quantum Harmonic Oscillator 
of Hamiltonian : 


- P ^ 02 A 

H g0) := + + Cp ^ 

when compared with the quantum evolution induced by Hpt) (g = 0), in the 
case where / is a T-periodic function and g a real constant. The reference 
(initial) state is taken to be an arbitrary “generalized coherent state” in the 
sense of Perelomov. We show that, starting with a quadratic decrease in time in 
the neighborhood of t = 0, this quantum fidelity may recur to its initial value 
1 at an infinite sequence of times . We discuss the result when the classical 
motion induced by Hamiltonian H 0 (t) is assumed to be stable versus unstable. 
A beautiful relationship between the quantum and the classical fidelity is also 
demonstrated. 
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1 Introduction 


A growing interest has been devoted recently on the study of the so-called “Quan¬ 
tum Fidelity” 

F g (t) := (Uo(t,0)ip,U g (t,0)ip) (2) 

for some reference state ijj that we take as initial wavepacket, Uo(t, 0) being the 
quantum unitary evolution operator induced by some Hamiltonian , and U g (t, 0) the 
quantum evolution for a perturbation of it, g being the size of the perturbation. The 
long-time behavior of F g (t) is of particular interest, and it has been studied for a 
large class of (time independent) Hamiltonians in a more or less heuristic way, and 
this long-time behavior has been suggested to depend sensitively on the regular 
versus chaotic motion of the underlying classical motion (see the bibliography). 

In most of these heuristic works (ID, i, i, ®, 0, HD, H2I, HD, HD, HD, 

HD, H2, HD, HD, HD, H3, HD, HD, E3, Ha, HD, HD, HU, HD, E9, E3, Ha, 

H21, gT], HU, [IS])? if is claimed that the Quantum Fidelity decays very fast to zero 
as time grows, when the underlying classical (unperturbed) dynamics is generically 
chaotic. 

Although the short time decay of the fidelity is rather well understood (33]), the 
arguments put forward in the above cited works are not entirely convincing, since 
they are either purely numerical, or extrapolate the “short t im e” behavior to guess 
the (Gaussian or exponential) decay at oo. 

In other approaches (|Z5|, JtjJ, (12]), the case of integrable or regular systems is 
considered as well, and seems to indicate the occurrence of an anomalous power law 
decay. Moreover, in a case of nearly integrable system a recurrence of fidelity has 
been exhibited; in this case the Quantum Fidelity manifests recurrences very close 
to the initial value 1 as time evolves. 

Thus it is a very intriguing subject of high interest to know more about the 
complete time behavior of the Quantum Fidelity when the underlying dynamics is 
chaotic versus regular. To our knowledge, no rigorous approaches of this topics have 
been attempted yet. It would be highly desirable to have a mathematically explicit 
description of the long time behavior of the Quantum Fidelity, although it is a very 
difficult task. 

This paper is a first attempt towards this rigorous understanding. It relies on 
a rather simplistic class of Hamiltonians for which the perturbed as well as unper¬ 
turbed quantum dynamics can be explicitly solved in terms of the classical dynamics. 
Moreover the reference quantum states are taken in a suitable large class of quan¬ 
tum states known as “generalized Coherent States” in the sense of Perelomov. The 
results are rather surprising : 
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- for a large class of reference states, the Quantum Fidelity never decreases to 0, 
but instead remains bounded from below by some constant. 

-in the unstable case, the Quantum Fidelity either decays exponentially fast to some 
non-zero constant, as t —> oo, or manifest strong recurrences to 1. 

-in the stable case the Quantum Fidelity always manifests strong recurrences to 1 
as time evolves. 

These facts are strongly related to the underlying SU( 1,1) structure underlying 
the corresponding Hamiltonians, as can be seen from the work of Perelomov ([23]). 

We also are able to show a strong relationship between “Quantum and Classical 
Fidelity” for this specific situation. 

To complete this Introduction, let us notice that a notion of Classical Fidelity 
that “mimics” the Quantum Fidelity has been proposed in the literature (El, 0, 
EH), where decay properties similar to those of the Quantum Fidelity appear. Thus 
it would be desirable to understand more deeply the relationships between the Clas¬ 
sical and Quantum Fidelity on a firm mathematical basis. We shall pursue this 
investigation in future publications, notably in the semiclassical limit (see Ena, 0). 


2 Calculus of the Quantum Fidelity 

Let us consider the following operators in 7i = L 2 (M) with suitable domains (see 


Kn = 


Q 2 + P 2 g- 


2 Q 2 


r - Q 2 -P 2 .Q.P + P.Q g 2 
K± = -:- Ti- 


4 ' " 4 2 Q 2 

where Q is the usual multiplication operator by x and 

P d 

p := -'Tx 

K 0 and K± satisfy the usual commutation rules of SU(1, 1) algebra, namely 

[K 0 ,K±]=±K± 


(3) 

(4) 


(5) 


[K_,K + ] = 2K 0 (6) 

We define the “generalized coherent states” (squeezed states) as follows : given 


any (3 G C 


Sg = exp (j3K + — 0K_) 


(7) 
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0/3 = S(P)ll) 

0 being a normalized state in 7i such that 

Ji_0 = 0 


( 8 ) 

(9) 
( 10 ) 

2 ■ v 4 " J (U) 

We shall focus on the following cases g = 1, g — \/3, where 0 has the following 
form : 


2/t o 0 = (a + -)0 


with 


a := ~ + t/y + 2# 2 


with 


</j(x) = cix 2 e x2/l2 

(12) 

x{x) = c 2 x 3 e ~ x2 ^ 2 

(13) 


(14) 

-i/I 

(15) 


(omitting the factors in tt). 

This makes ip and x to be (normalized to 1) finite linear combinations of Hermite 
functions. 

It has been shown (see |24j ) that for 0 = ip, 0^ has the following form : 

0/3 = cnrV 2(u ~ e) exp i((re~ ( “~ e) ) 2 ^ (16) 

where the constants u, it, 9, e are suitably determined from (3 € C , whereas for 

0 = X, "0/3 is : 


= c 2 x 3 e exp ( ~ + 0^ - hu - e) - i 


(17) 


Now the evolutions of 0^ with respect to U\(t, 0) and U^(t, 0) respectively, to¬ 
gether with 0) are completely explicit. 


4 





Consider a complex solution of the classical equations of motion induced by 
Hamiltonian Ho(t) : 

x + fx — 0 (18) 

and look for it in the form 

x = e u+ie (19) 

where the functions 1 1 —>• u and 1 1 —>• 0 are real. 


We assume for u(t) and 0(f) the following initial data : 

u( 0) = u 0 
h(0) = M 0 

0 ( 0 ) = 00 

0(0) = e - 2 (“o-0 

Since / is real, the wronskian of x anf x is constant and equals 

2 ie 2e = 2 i6(t)e 2u{t) (20) 

This yields 

x = [u + id T (h T i0) 2 ]x = [h + m 2 — e~^ u ~^]x = —fx (21) 

and therefore u obeys the following differential equation 

m + m 2 - e“ 4( “" e) + / = 0 (22) 

We have the following result : 


Lemma 2.1 Let H g = p fp- + p. Then the quantum propagator for H g (t ) is of 
the following form : 


U g (t 0) = e iil Q 2 / 2 e -i(u-£)(Q-P+PQ)/‘ 2 e -i(e-0o)H ge i(u 0 -e)(Q.P+P.Q)/2 e -iu 0 Q 2 /2 ( 23 ) 

The same formula holds for the propagator f/ 0 (£, 0) of Ho(t) with H g replaced by Hq. 
Proof : Let us denote 


V(t) ’= /2 e -i(u-e)(Q.P+P.Q)/2 e ~i6H g 


(24) 


We have : 


ijYO) = + |e“ q2/ \Q.P + P.Q)e- luQ2 ^ V(t) 


(25) 


5 



( 26 ) 


^Q e iuQ 2 /2 e i(u-e)(Q.P+P.Q)/2jj^ e -i(u-e)(Q.P+P.Q)/2 e -iuQ 2 /2' s j y(^) 

and since 

e iuQ 2 /2p e -iuQ 2 /2 = p _ 

we have : 

e™ Q2/2 (P.Q + Q.P)e~ iilQ2/2 = P.Q + Q.P - 2 iiiQ 2 
Therefore the hrst line in the RHS of (25) is 

(~-u 2 )Q 2 +^(Q.P+P.Q)^V(t) 

Furthermore 


(27) 


(28) 


e -i{u-e)(Q.P+P.Q)/2p e ~i{u-e){P.Q+Q.P)/2 _ f JLS\ e 2 ( u ~ e ) _|_ Q_ e ~ 2 ( u ~t) ( 29 ) 

\ 2 Q / 2 


which implies that the second line in the RHS of (25) is 

where we have used that 


V(t) 


(30) 


6 = e - 2{u ~ e) 

Collecting the different terms, we get : 

. d 


S VW =(t + S + T + (31) 


' p 2 a 2 O 2 ' 

= ^ T + Q2 


using (22). □ 


We shall now consider the “quantum fidelity” in two particular cases g = 1 and 
g = -\/3, starting respectively with the initial states ^i = S(/3)ip , and = S(/3)x ■ 

Fi(t) = (U 0 (t, 0)^/3, 1 , C s (t, 0)^/3, 1 ) (32) 

and similarly for F 2 (t) with t/jpp replaced by 7 / 53 , 2 - 
Theorem 2.2 (%) Let 5 = 1 and 7 /; ,3,1 = S(/3)y. Then we have 

Fl (t) = l + l e - 2 mt)-m) (33) 

3 3 

(ii) Let g = and ipp j2 = S(P)x- Then we have : 

F 2 (t) = - + (34) 

5 5 
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Proof : Since x 2 is expanded as 


x = ^H 2 (x ) + -H 0 (x) 


it is clear that 


* .2 1 


Uo(t, 0)^/3 ,1 = Cl exp - ~(u(t) - e) 

X |i e -SW ( t )/2 ^ 2 ( a;e - ( «(t)-e ) ) + I e -^(t)/2-2i0(O) | exp 


x 2 e -2(u(i)-e) 


and 


U g {t)il)p, i = V(t)<p 

= Cl exp ( -^i0(f) + - I(a;e- (u(t) - e) ) 2 


from which we deduce that 


( 35 ) 


(36) 


Flit) = c\ 

4/1 1 


1 p -2i(e(t)-«(o)) 

x\x 2 --) + ---x 2 


e~ x dx 


+ _ e -2im-m) ) = _ + i e -2*(e(t)-e(o)) 


3 V2 4 / 3 3 

(ii) follows from a very similar calculation using that 


(37) 


H-dx) = 8x 6 — 12x 


□ 


Let us now assume that g has an arbitrary real value (not specifically of the form 

g = \fk[k + l)/2 for k G N which gives rise to integer values of a — \ + \J\ + 2 g 2 ). 

Then the state 0 := cx a e~ x2 t 2 is no longer a finite linear combination of Hcrmite 
functions. It has nevertheless an infinite expansion on the basis of Hcrmite functions 

^ Pn • 


oo 

0=A n (/? n 

n=0 

and one can establish the following general result about the corresponding “quantum 
fidelity” : 

Theorem 2.3 

OO 

F g (t) = exp (~ia(6(t) - 0(0))) |A n | 2 exp(m(0(f) - 0(0))) 

n =0 
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which, since Yl^=o |A n | 2 = 1 , returns in absolute value to 1 as long as Off) —0(0) = 
0 (mod 2n). If a = p/q is rational, then the quantum fidelity recurs exactly to 1 (not 
only in absolute value) provided that Off) — 0(0) = 0 (mod 2q7r). 

Proof : Again according to ref. [21] . we have 

'fig := S((3)fi> = e - i( - a+1 i )e °D(u 0 )fi 
where by D(u) we denote the following unitary operator 


Then 


D(u) • — g»«Q 2 /2 i e ~i(u-e)(Q.P+P.Q)/2 


whereas 


U Q (t,0)i/>p = D(u t )e- l ^- eo ^ 0 J2^e- ieo{a+ ^^n 

n =0 


= D(ut) A ne -<fi*(n+ 3 )-H fl b(»-«) ¥ ,, 

72—0 


U g (t,0)i/;p = e- ie *^D(u t )^ 

so that : 


(Uo(t, 0)lpp, Ug(t, 0)fi>p1pp) = AnA me ^(n+i)-^o(n-a)-^( a+ I)(^ ) ^ 

n,m 


□ 


oo 


= En 

72—0 


2 ^i(6t—do)(n—a) 


3 Discussion of the Result 

Since / is a T-periodic function, Floquet analysis applies to equation (18) which 
is nothing but the well-known Hill’s equation. Depending to the parameters char¬ 
acterizing the function / (as for example 7 and 6 in the case of Mathieu equation 
where f(t) =7 + 6 cos cut), the solutions can be either stable or unstable. In all 
cases the quantum fidelity iq ft) and F 2 (t) are bounded from below by some positive 
constant in absolute value, and therefore never decrease to zero. 

The phase Off) is determined by 

Off) - 0(0) = e 2e [ e~ Ms) ds 


0 


(38) 


In the stable case, the function t i—► u(t) is T-periodic. Therefore &(t) grows 
from — oo to +oo as time evolves from — oo to +oo. It follows that there exists an 
infinite sequence of times {£*.} where F\ (t) recurs to its initial value 1, and an infinite 
squence {t' k } where \F\ (t) attains its minimum 1/3. The same statement holds for 
F 2 (t) where 1/3 is replaced by 1/5. 



In the unstable case there is some positive Lyapunov exponent A and some real 
solution of Hill’s equation such that 

x(T) = e XT x(0) 
x(T) = e XT x(0) 

Here, since we deal with complex solutions x(t) such that the Wronskian of x 
and x is non zero, we deduce that |x(t)| > 0. Moreover, depending on the instability 
zone, \x(t)\~ 2 = e _2 “^ can be integrable near ±oo , or diverge at ±oo. These topics 
are clearly detailed in ref. m- In the first case the conclusion is that there exists 
two constants 6± such that 


6{t) —> 9± as t —> ±oo (39) 

This is the case for the Inverted Harmonic Oscillator (/ = — 1) that we describe 
in the last Section. 

Therefore the Quantum Fidelity in this case behaves generically as described in the 
following picture : 
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1/3 


t 

If the instability zone in which the solution xft) lies is such that |x(t)Q 2 is not 
integrable near ±oo, then the figure is typically similar to that of the stable case, 
which shows infinite recurrences to 1 of the Quantum Fidelity. 

4 Link with the “Classical Fidelity” 

We can call “classical infidelity” the discrepancy between two classical trajecto¬ 
ries, along their evolution, when they merge from the same initial phase-space 
point at t — 0. The classical fidelity is thus the possible crossing of the two trajec¬ 
tories governed by H 0 (t) and H g (t), in phase-space, as time evolves. 

From now on we assume that e of Section 2 is determined by (g > 0) : 

e 2e = gV2 (40) 

Proposition 4.1 Let x(t) be a real trajectory for the time-periodic Hamiltonian 
H 0 (t), and y(t) a real trajectory for Hamiltonian H g (t). Assume in addition that 
they satisfy x(0) = y(0) G M, ±(0) = 2/(0) e M. Then we have : 

I x(t) - yft) | = | y(t) | (1 - cos 0(t)) (41) 

where Oft) is defined as follows : 

9{t) \= gy/2 j y(s)~ 2 ds = Oft) — 0(0) (42) 

Proof : Let zft) be a complex solution of equation 

z + fz = 0 (43) 

It can be written as zft) = e u+ld as above, where t h u and t e-> 0 are real 
functions. Assume that xft) is such that x(0) > 0. Then it is easy to see that 


10 







x(t) = 3f lz(t) = e u ^ cos 0(f) is a solution of equation (43), and that the positive 
function y(t) := e u ^ is a solution of 

2q 2 

y + fy — 3- = 0 (44) 

ir 

which means that it defines a classical trajectory for H g (t). 

Namely we get from equ. (42) that 



^-9{t) = gV2e~ 2u(t) 

CLL 

(45) 

and therefore 



x = (it + u 2 — 2 g 2 e 4u )x = —fx 

(46) 

so that 

(u + u 2 - 2g 2 e~ 4u + f)e u = 0 

(47) 

which is nothing but equ. (44) noting that 



V = (u + u 2 )y 

(48) 

Moreover, they satisfy 

®(0)=J/(0), ±(0) = 2/(0), and 



|®(f) - y(t) 1 = y(t)( 1 - cos 0(f)) 



In the case where x(0) < 0, we just take x(t) = —ffiz(t) and y(t) = —e u ^\ which 
completes the result. □ 

Conclusion : The classical infidelity vanishes for 0(f) = 2kn (k G Z), which 
precisely gives rise to recurrences to 1 of the “quantum hdclity” (Theorem 2.3). We 
expect this remarkable property to be true in more general situations, in particular 
in the Semiclassical Regime (see 

By “vanishing of the classical infidelity” (thus “classical hdclity”) we meant 
that given any solution of the unperturbed dynamics, there exists a solution of the 
perturbed one that coincides with it at the origin, and at any values of time solving 
the equation 


gV 2 / ds y(s ) 2 = 0 (mod 27r) 


One may ask whether this holds true for general solution x(t), y(t ) of equ. 
(43-44). This is answered in the following particular cases below. 

• Particular case / = 1 
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(The case / = u 2 could be treated as well.) 

When we apply the result above, we find that solutions x(t) = A cost and y = A 
coincide for t = 2kir, fceZ, and that y obeys equation (44) provided A = ( 2 g 2 ) 1 / 4 . 
One shows that a more general result holds, involving general solutions of the 
equations under consideration. 

The general solution of equ. (44) is of the form (apart from the sign before the 
square root) : 


y{t) = v/aT (3 cos 2t + 7 sin 2t (49) 

with a, (3, 7 related to each other by the relation 

« 2 - (£ 2 + 7 2 ) = 2 g 2 (50) 

It obeys the initial data 


2/(0) = \/a + (3 


m = 


7 


\J a + /3 


and the conserved energy is simply a. 

We take as real solution of equ. (43) (Harmonic Oscillator) : 


x 


(' t) = \] a + (3 cos t + 


7 


: sin t 


(51) 


\J OL + (3 

Which has the same initial data as y(t). Both functions being 2tt- periodic, the 
generic “classical infidelity” vanishes when t = 2/c7t, fceZ. 


• Particular case / = — 1 

(The case / = —cu 2 could be trated as well) 

Any complex solution of the differential equation (Inverted Harmonic Oscillator) 

z-z = 0 (52) 


can be written in the form 

z(t) = (a + ifye 1 + (c + id)e“* 
where a, b , c, d are real constants. We define : 

Z(t) :=\z(t )| 2 


( 53 ) 
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( 54 ) 


One can prove that y(t) = \JZ(t) obeys the differential equation 


2 g 2 n 
y-y —3- = o 

yi 


with 


g 2 = 2 (ad — be) 2 > 0 

It is important to note that this implies Z(t) > 0, Vt. 
Define x(t) := y{t) cos 9(t), with 0(0) = 0. Then, clearly 

x(0) = 2/(0) 

i(0) = y( 0) 

We want x(t) to be a real solution of equ.(52). Since 


(55) 


x — x = 


(z_P 


y/Z \ 2 4Z 
and since Z obeys : 


z - o 2 z\ cos e 


Z Z 2 „ 2 g 2 

- Z -— = 0 

2 4Z Z 


9 Z \ 

—j= + Q\/~Z j sin 6 (56) 


(57) 


the RHS of equ. (56) vanishes provided 


0 2 = 2 g 2 Z~ 2 

Then x(t) and y(t) are two trajectories for Hamiltonians p and 
respectively, with the same initial data, and we have 

\ x (t) -y(t )I = y(t)( 1 - cos $(t)) 

Using the particular solution of equ. (52) 


(58) 


P 2 -Q 2 , £_ 
2 + Q 2 


(59) 


z[t) = 


1 + i t 1 - i 


that can be rewritten as : 


z(t) = V cosh 2 1 exp(- arcsin tanh 2 1) 

i.e. with u(t) = 4 log cosh 2t (which satishes u(0) = h(0) = 0), the formula (23) 
can be rewritten in the simple form 


U g (t, 0) = exp ( -Q 2 tanh2tJ exp ^-(Q.P + P.Q) log cosh 2tj exp arcsin tanh 2t 

(60) 
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which holds true for g — 0, and g = l/\/2. For g = 0, this is nothing but the 
well-known Mehlcr’s Formula. 

Study of 9{t) 

Recall that the fidelity F g (t) strongly depends on the reference state via the 
constants a, 6, c, d 


m = 



ds 


o (a 2 + b 2 )e 2s + (c 2 + d 2 )e~ 2s + 2 (ac + bd ) 


(61) 


= ( arctan 


ac + bd + (a 2 + b 2 )e 2t 
| ad — be | 


arctan 


ac + bd + (a 2 + b 2 ) 
| ad — 6c | 


This implies that d(f) —> 9± as t —> ±oo exponentially fast in the future and in 
the past. The calculus is especially simple in the particular case where we choose 
a = d = g/V%, c = b = 0 : 


9(t) = arctan(e 2t ) — — (62) 

No classical fidelity happens, and the square of the quantum fidelity in the case 
g — 1 is nothing but 


as t 


±oo. 


m i 2 


5 8 

9 9(e 2t + e~ 2t ) 




9 9 


(63) 


Here the symmetry is perfect between the future and the past, and the Quantum 
Fidelity has no recurrences to 1, but instead decays exponentially fast to 5/9. 



5 CONCLUSION 

It has been suggested in many recent works that the type of decay of the Quan¬ 
tum Fidelity may help to discriminate between chaotic or regular underlying clas¬ 
sical motion; in other words the hypersensitivity of quantum dynamics under small 
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perturbations, as measured by the type of decrease of the Quantum Fidelity, could 
be a signature of what is often called (rather improperly) “Quantum Chaos”. 

Thus it is highly desirable to have a better understanding of how this function 
of time (represented by equation (2) for suitable class of reference quantum states 
0 ) behaves at infinity in general as well in particular systems. In this paper we have 
been able to describe the full time-behavior of the Quantum Fidelity for a rather 
specific class of systems, and for reference states in a suitable class of “generalized 
coherent states. The underlying SU(1, 1) structure possessed by these systems allows 
us to perform an exact calculus of the Quantum Fidelity, and to compare it with 
the “Classical Fidelity” of the corresponding classical motion. This classical motion 
can be either stable, or unstable with positive Lyapunov exponent. The Quantum 
Fidelity has the following remarkable properties : 

-either it decays to a (generally non-zero) limit in the past and in the future 
-or it manifests an infinite sequence of recurrences to 1 as time evolves. 

This sheds a new light on this question which has been addressed in a great 
variety of cases in the physics literature, and where the Quantum Fidelity is generally 
claimed to decay very rapidly to zero. Thus the first mathematical study presented 
here on the long time behavior of the Quantum Fidelity could allow in the future a 
better understanding of these features in more general situations on a rigorous level. 
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important information about solutions of Hill’s equation (see ref 0). 
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